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^ Abstract 

Analyzing the dilepton production from out of equilibrium quark-gluon 
plasma, Le Bellac and Mabilat have recently pointed out that, in the 
reaction rate, the cancellation of mass (collinear) singularities takes place 
only in physical gauges, and not in covariant gauges. They then have 
estimated the contribution involving pinching singularities. After giving 
a general argument for the gauge independence of the production rate, we 
explicitly confirm the gauge independence of the mass-singular part. The 
contribution involving pinching singularities develops mass singularities, 
which is also gauge dependent. This "additional" contribution to the 
singular part is responsible for the gauge independence of the "total" 
singular part. We give a sufficient condition, under which cancellation of 
mass singularities takes place. 
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In the past years, much effort has been made to incorporate quantum field theory with 
nonequilibrium statistical mechanics, among which, we quote those of Altherr and 
Seibert |J, Altherr p|, Baier et al 0, Le Bellac and Mabilat Q and the present author 
p. Out of equilibrium, pinching singularities appear |J in association with self- 



to the rate of hard-photon production from nonequilibrium quark-gluon plasmas is 
made in 0. A renormalization scheme of number density is introduced in ||, such 
that the pinching singularities disappear. Le Bellac and Mabilat Q are the first who 
have explicitly analyzed the infrared and mass (collinear) singularities in "lepton-pair" 
production rate. 

First of all, let us summarize the results of Q. The production rate of a lepton 
pair from a quark-gluon plasma is proportional to IT(Q) = —iHi 2 (Q), where Hi 2 (Q) 
is the (1, 2)-component of the photon self-energy part in real-time massless QCD. 
Following ||, we deal with H(Q) of a scalar "photon." To two-loop order, H(Q) 
receives two contributions. The one lis comes from the diagram with self-energy 
inserted quark propagator and the one LV comes from the diagram with photon- 
quark vertex correction, 



energy inserted propagator. It has been shown in |2| that resummation of self-energy 
part eliminates the pinching singularity (see also ||). An application of this result 



n(Q) 




(i) 



(2) 



3,1=1 




xS 2l (P - Q)rS n (P -Q-K)} Ay(K) , 



,feau ge) W £ ( _ ) i+/ Tr [Sij{P _ K )YS J2 (P) 



(3) 



where 



^l(P) = ^9\-y +l l^9^ use \K)rS jl (P-K)YA, l (K). 



(4) 



The (part of the) gluon propagator Ay (if) takes the form Q 



A u (if) 
A 12 (K) 
A 21 (K) 



A* 2 (K) = (1 + f(K))A R (K) + f(K)A A (K) , 

f(K)(A R (K) + A A (K)), 

(l + f(K))(A R (K) + A A (K)), 



(5) 
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where A^a(K) = ±i/(K 2 ±iO + e(ko)). The quark propagator Sij(K) takes the form 
Sij(K) = l£Aij(K), where Ay (IT) is given by fl5|) with the substitution / — > —f. f 
(f) is related to the distribution function of gluon n (quark ft) through 

f(K) = -6(-k ) + e(k )n(\k \,e(k )k), 

f(K) = e(-k )+e(k )n(\k \,e(k )k), (6) 

where k = k/fc with k — |k|. Form of g^ nge \K) in (^) and (f|) depends on the gauge 
choice: 

g^\K) = g^- v ^^ (covariant gauge) , (7) 
g$(K) = g^-h( K ^ + n^K„) + 1 ^ (Coulomb gauge) , (8) 
where n M = (1,0). Observe that Q 

X; 5 1J (P)S, 7 (P)^ 2 (P) = + P^) , (9) 

0,1=1 

F (n) (p) = _/(p)( A 2 (p) _ Ai(P))^ReSn(P)^ 

- ±/»(A* (P) + A^(P))/(S 12 (P) - £ 21 (P))^ , (10) 
F«(P) = A R (P)A A (P)/[(1 - /(P))£ 12 (P) + /(P)£ 21 (P)]/ . (11) 

p( ra ) is the "normal term," which is the counterpart of the one that is present in 
equilibrium thermal field theory (ETFT), while is the "pinch term," which is 
absent in ETFT. Substituting ([5|) - ( fTTf) into (0), we have, with obvious notation, 

n s (Q) = n^(Q) + n^(Q). 

Le Bellac and Mabilat |4| have shown that, in Coulomb gauge, the mass singularities 
cancel out both in (Q) and Ily(Q), Eq. (^|). While, in the covariant gauge, the 
cancellation holds only for IIj? (Q), and in Hy(Q) there survives mass singularity. 
Then, the authors of have concluded that whether or not the singularity cancella- 
tion takes place is gauge dependent. [In Appendix A, we show in a gauge- independent 
manner how the cancellations of mass singularities take place in IT^ 1 (Q), reconfirming 
the result of 0.] 
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We first verify that H(Q), Eq. ([!]), is gauge independent. Proof goes just as in 
vacuum (T = 0) theory. Consider the difference 



sn(Q) = n E (Q) -n s (Q) , (12) 

covariant Cou 

where "Cou" stands for Coulomb. Observe that, from ([/]) and ([8]), the difference 
g^° v \K) — g®(K) is proportional to and/or K v . Then, in evaluating 5U(Q), we 
can use Ward-Takahashi relation, 

S jk (P - K)tfS u {P) = -i(-) h 5 k iS jk (P -K) + i(-) j 5 jk S kl (P) , 

with no summation over repeated indices. After doing this, we see that, among many 
terms in the resultant expression for SH(Q), complete cancellations occur, so that 
SU(Q) vanishes and then is, of course, free from mass singularities. 

On the light of the above observation, let us make a closer inspection of the 
results of ||. As a covariant gauge, as in ||, we take the Feynman gauge (r) = in 
(0)) throughout in the sequel. We analyze IT S (Q), Eq. (Q), which includes Q^ v = 
f^^f — Simple algebra yields 

g^gr = -2P 2 g, (13) 
6g, u g^ = 2P 2 ^ 2 [(2p -k )^-p Q g]+O((P 2 ) 2 ), (14) 

where Sg^ = g^ v — gj®(K) and use has been made of K 2 = (P — K) 2 = 0. Thus, 
both g^vQ^ and bg^Q^ are proportional to P 2 . Now we notice that P 2 Ar (P) Aa(-P) 
= P/P 2 . This means that the pinching singularity in the "pinch" term FW, Eq. ([TT|) , 
turns out to be a mass singularity. The (P 2 ) 2 term in (|14D does not lead to mass- 
singular contribution. As in |4[], let us restrict our concern to singular contributions 
and ignore the (P 2 ) 2 term. 

As far as mass-singular contributions are concerned, above observation 5U(Q) = 
tells us that II{> mg | Fc = — , where "Fey" stands for Feynman. ny mg | Fe with 
Q = (<?o? q = 0) is explicitly evaluated in ||. As a check, in Appendix B, we evaluate 
<m^ sing (g , 0) and confirmQ 5U^ sing = -n^ ing | Fey . 



'''There is a missing term in IIv|p m @] (see Appendix B). 
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The singular part of n(g = 2k, 0), being gauge independent, is also evaluated in 
Appendix B, 



■QSing 
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Cou 
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+ n (p),sm 5 
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/■°° z 2 + 2 f 
+ / az n{KZ, k) 

JO -2 



o 



ct2; F — : — n{Kz } k) 



16 2l 1 rdQs„ 
H aa s K In — / — -n(K, — k) 

37T 



z 2 - 1 

1 ( z _ 1)2 _|_ i 

dz ti(kz, k)n(/t(l — z), k) 



+ / dz- i 7i(«3,k)(l - l),k) 

Ji z 

poo ( z _|_ X) 2 + 1 

+ / dz- (1 +n(/cz,k))n(«(l + z),k) 



(15) 



Here cif^ is the element of the solid angle in the k-space. The cutoff factor e y is 
defined by 1 — |p • k| > e y (cf. Eq. ( |B.1| ) in Appendix A) and e z is the infrared cutoff 
k > e z K. 

Let us clarify the relation between the present result and the result of @]. We 
start with picking out from F^(P), 



E(P) = (1 - /(P))E 12 (P) + /(P)E 21 (P) . 



(16) 



For the purpose of estimating , Le Bellac and Mabilat Q have analyzed E(P) 
within the hard-thermal-loop resummation scheme. The net production rate of an 
(anti) quark is given by T(P) = Tr[iS(P)P]/(4p), with P = (p, p) for quark and 
P = (— p, — p) for antiquark. Arguing that T(P) on the mass shell P 2 = 0, being 
gauge independent, is relevant to I1 S (Q), the authors of || have concluded that n(Q) 
is gauge dependent since n sing (Q) is. It is clear from the above argument that this 
is not the case. As has been discussed above in conjunction with (|14T ), i^(5S(P)P 
(as well as P'S(P)/ > | F ) vanishes on the mass shell P 2 = 0. [As a matter of fact 
£12(21) (P) in ( |16D vanishes on the mass shell P 2 = 0, since P 2 , K 2 and (P — K) 2 
cannot vanish simultaneously] However, as noted above, in calculating 8F^ (cf. Eq. 
(0)), the factor P 2 , Eq. (0), "eliminates" one A and the mass-singular contribution 



6U 



(p) > sing 



fas well as n 



(p), sing I 



I Fey 



) emerges. Thus we have learned that the mass- 
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singular contribution does not come from the gauge-independent quantity, 



^(^| p2 = = - 2 «^| p2 = (= °) 

but comes from the gauge-dependent quantity 

d^(P)fi/dP 2 \ = PH(P)P/P 2 



\p 



The observation made above in conjunction with (14) applies to the contribution 
(to <5F (p) ) from the soft-X region (cf. fly]) with (D), in which K 2 ^ 0. [The soft 



(P— i^)-region is not important, at least, for the system, which is not far from thermal 
and chemical equilibrium.] Above observation on (ftl^ sms holds as it is, except that 
Si2(2i)(-P) does not vanish on the mass shell. It is to be noted in passing that n^| Fe 
develops pinch singularity. This is because, in the present case, K 2 ^ 0, we have, in 
place of (|l|), g^Q^ = -1P 2 H. + 2K 2 f. The second term on the right-hand side 
leads to pinching singularity in rij? . Because of the factor K 2 , which is small, the 
"residue" of the pinching contribution is relatively small. 

Finally, for an instructive purpose, we derive a sufficient condition for the absence 
of mass singularities, i.e., the condition under which flT5| ) vanishes. The condition is 

[1 + f(K)][l - f(P - K)] = [1 - /(P)][l + f(K) - f(P - K)] , (17) 

which holds for equilibrium case. For simplicity of presentation, we take the self- 
interacting scalar theory counterpart ||] to (|T7|), 



[1 + f(K)][l + f(P - K)} = [1 + /(P)][l + f(K) + f(P - K)} . (18) 

The same conclusion, obtained below, applies to ([17]). 

We start with the following observation. It is well known that the reaction rates 
are free from mass singularities if all the (energy) degenerate states are prepared in 
the initial state (see e.g., [0])). Mass singularities arise from collinear configurations of 
massless particles. For convenience, we enclose the system into a box of volume V (= 
L 3 ) and employ a periodic boundary condition. In the interaction representation, the 
scalar field may be expressed as a superposition of plane waves, whose "coefficients" 
a^. and at are the annihilation and creation operators, respectively. Here k,- is the 
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momenta of a particle with mode j. The, above observation tells us that for the 
system characterized by the density operator 



be remarked that, even in the present context, the form of the density operator takes 
much more complicated form, since the system is in general not uniform in space- 
time. However, as far as a reaction rate taking place in the system is concerned, 
all the propagators describing the reaction rate have @, |6|] common center-of-mass 
or macroscopic coordinates of the space-time region, where the reaction takes place. 
Then it is sufficient to use (|19"1). 

Let us show that (|T§| ) leads to (|TJ). We start with computing Trp 




(19) 



reaction rates are free from mass singularities. In (|T9|), kj = |kj] is the energy of a 
particle with kj and /(kj) is a (smooth) function. If /(kj) = 1, the argument of p 
in (|T9l ) is essentially an energy operator in the initial (remote past) state. It should 



Trp 





where 



i 




Now, the number density is evaluated as 



n{k,ki) 




Armed with the above preliminary, we take (|T8|) with po > ko > 0. The right-hand 
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side of (|T8| ) reads* 

[1 + n(k n , k n )] [1 + n(kj, k,) + n(\k n - kj\, k n - k,-)] 

(20) 

where iV(£, k) = l/^e 1 ^ 9 ^ — 1). Since we are interested in collinear configuration, 
k n ~ k 3 ~ k n — kj and | k n — kj | ~ k n — kj . Then we have 

Eq. J^e^e^^N^k^N^K-k^g, 

which is the left-hand side of (0). For other regions of (po, ko) than the one studied 
above, one can similarly prove (|18|). 

It should be noted that, since p(E) and g(k) are arbitrary functions, the resultant 
form for n(k, k) "covers" a wide class of functions. 
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Appendix A Absence of mass singularity 

in nf >(Q) 

In this Appendix, we show that II^ 1 (Q) is free from mass singularities, reconfirming 
the result in ||. Manipulation goes as follows. Substitute fW, Eq. QTD|) , into (0) 
with (|). Use the form of gj^ uge) , Eq. © or Eq. (§, and forms for A Lj and S h (cf. 
Eq. (P)). The resultant expressions may be rearranged as 

•tFor the mode-overlapping "points", k„ = kj, k„ = 0, kj = 0, ( po[ ) does not hold in general. The 
contributions from such "points" are 0(1/ V) smaller than the "reference contribution" and vanish 
in the limit V — > oo. Incidentally, for the case of equilibrium system, (|20| ) holds B. 
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X( 



X 



SOw)/(p)(i _ /(p _ Q))5 e ((P - Q) 2 



- + /(Jf) 5 e (X 2 ) (A* (P)A fl (P - K) + A\(P)A A (P - K) 



5 - /> " iO) U(P - Kf) (A%(P)A R (K) + A\(P)A A (K 



(A.l) 



where £? Ml/ is as in fll3D and 5 e (A; 2 ) = e(k )5(K 2 ) etc. For the Coulomb gauge, 
^feauge) _ _(t)_ p Qr foe covar j an t gauge, <^f, auge - ) = g^ v + ^K^K^d/dK 2 ), where 
d/dK 2 applies to 5(K 2 ), Ar(K) and A A (K). Equation ( |A.1| ) is manifestly free from 
mass singularities. Mass singularity arises from the terms A 2 r , a AP)Aa{r)(P — K) and 



A|j/^(P)A^(^)(J^). In obtaining (|A.1|) , cancellations occur between those terms. 



It is also obvious from (|A.l ) that 11^ is free from divergence due to infrared 
singularities, provided that, as k — > + , f(k, k) oc A;~ n with n < 2 and /(&, k) oc A;~ n ' 
with n' < 2. In actual computation of II, for a propagator with soft momentum, one 
should use hard-thermal-loop resummed effective one. 



Appendix B Computation of the singular part of 

nf (Q) 

Here, we compute the singular part of the "pinch" contribution n|? (?o,q = 0). 
Substituting (|TTD into (|2|) with (|) and using (p~3f ) and flU]), and the forms for Ay and 
Sy, we have 



nr ,Smg = - l ^aa s Q 2 f^f ^Ld^\l - f(P - Q))5 e (K 2 )5 £ ((P - Kf 

D7T J 271 J 2TT 

1 



x6 e ((P - Q) 2 ) P T2 [f(K)f(P — K) — /(P)(l + f(K) - f(P - K))] , 



where 

G (Fcy) 
q(Cou) 



ko 

(go - kp) 2 + k 2 
2q k 
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Here use has been made of po = qo/2, which comes from (P - Q) 2 = 0. Making the 
change of variable P — > P + K, we extract the mass-singular part, 

= tt(5(P 2 ) fdkk 2 f dk 5{K 2 ) f Ey d(p ■ k) ' 



-i+e y poh - p ■ k 

-6(P 2 ) fdkk I dk e{ Po k )5{K 2 ) In — . (B.l) 
p J J e,. 



Using cutting off the infrared region, (= kz) > e z n, and changing the inte- 
gration variable suitably, we arrive at the final form. As it should be, ^II^ ms = 
n^' sins | — IIj? mg | is equal to — Ily ng , which has been evaluated in [ffl. There 

I Fey I Cou 

is a missing term though in ITy ms in [ffl, 



aa s K 2 (n(K)) 2 In — f dz 

STY dy J €y Z 



where k = qo/2. The form for n smg , being gauge independent, reads flT5| ) in the text. 
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